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We report the first measurement of the forward-backward asymmetry and the ratios of Wilson 
coefficients in B —* K*£ + £~, using 386 xfO 6 BB pairs that were collected on the T(4S) 
resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. We also 
present the first evidence of the decay B~ — » t~v t , using 4f4 fb -1 of data. 



1 Introduction 

Flavor-changing neutral current b — > s processes proceed via loop diagrams in the Standard 
Model (SM). The b — > s processes are sensitive to new physics effect. If new heavy particles can 
contribute to the decays, their amplitudes will interfere with the SM amplitudes and thereby 
modify the decay rate as well as decay distributions. Such contributions may change the Wilson 
coefficients 1 that parametrize the strength of the short distance interactions. The b — > sl + i~ 
amplitude is described by the effective Wilson coefficients C7 , Cg ff and (5fg- Measurement 
of the forward-backward asymmetry and differential decay rate as functions of q 2 and 9 for 
B — > K*£ + £~ constrains the relative signs and magnitudes of these coefficients 2,3 . Here, q 2 is 
the squared invariant mass of the dilepton system, and 6 is the angle between the momenta of 
the negative (positive) lepton and the B (B) meson in the dilepton rest frame. 

The purely leptonic decay B — > t~v t proceeds via annihilation of b and u quarks to a 
W~ boson in the SM. It provides a direct determination of the product of the B meson decay 
constant fs and the magnitude of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
\V u b\. The branching fraction is given by 
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where Gf is the Fermi coupling constant, ms and m T are the B and r masses, respectively, and 
tb is the B~ lifetime 4 . Purely leptonic B decays have not been observed in past experiments. 
The most stringent upper limit on B~ — > t~v t comes from the BaBar experiment: B(B~ — > 
r-p r ) < 2.6 x 10" 4 (90% C.L.) 5 . 

The Belle detector is a large-solid-angle magnetic spectrometer consisting of a silicon ver- 
tex detector, a 50-layer central drift chamber (CDC), a system of aerogel threshold Cerenkov 
counters (ACC), time-of- flight scintillation counters (TOF), and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that provides 
a 1.5 T magnetic field. An iron flux-return located outside of the coil is instrumented to identify 
K L and muons. The detector is described in detail elsewhere 6 . 

2 Measurement of Forward-Backward Asymmetry in B — ► K*£ + £~ 

We use a 357 fb _1 data sample containing 386 x 10 6 BB pairs taken at the T(4S) resonance. We 
also study the B + — > K + i + i~ mode, which is expected to have very small forward-backward 
asymmetry even in the existence of new physics 7 . The following final states are used to re- 
construct B candidates: K*°£+£-, K* + £ + £~, and K+£ + £~, with subdecays K*° -> K + ir~, 
K*+ -» K%tt+ and K + ir°, K° s -> tt+ti-, and tt° -» 77. Hereafter, K*°£ + £~ and K*+e+£- are 
combined and called K*£ + £~ . We use two variables defined in the center-of-mass (CM) frame 

to select B candidates: the beam-energy constrained mass M^ c = E 2 cam — p B and the energy 
difference AE = Eb — ^beam, where ps and Eb are the measured CM momentum and energy 
of the B candidate, and -Ebeam is the CM beam energy. The dominant background consists of 
BB events where both B mesons decay semileptonically. We suppress this using missing energy 
and cos9 B , where 6* B is the angle between the flight direction of the B meson and the beam 
axis in the CM frame. Backgrounds from B — > J/i/jX s ,iJj(2S)X s decays are rejected using the 
dilepton invariant mass. The signal box is defined as \M^ C — ms] < 8 MeV/c 2 for both lepton 
modes and —55 (—35) MeV < AE < 35MeV for the electron (muon) mode. 

We perform an unbinned maximum-likelihood fit to the M^ c distribution to determine the 
signal yield. The fit function includes signal, cross-feeds and other background components. 
In the fit, all background fractions except the dilepton background are fixed while the signal 
fraction is allowed to float. We obtain 113.6 ± 13.0 and 96.0 ± 12.0 signal events for K*£ + £~ 
and K + £ + £~, respectively. Figure 1 shows the fit result. 




Figure 1: Mbc distributions for (a) B — > K*£ + £ and (b) B — » K + £ + £ samples. The solid and dashed curves 
are the fit results for the total and background contributions. 

We use B — > K*£ + £~ candidates in the signal box to measure the normalized double differ- 
ential decay width. For the evaluation of the Wilson coefficients, the NNLO Wilson coefficients 
Ci of Ref. 8 are used. Since the full NNLO calculation only exists for q 2 /ml < 0.25 region, we 
adopt the so-called partial NNLO calculation 9 for q 2 jm\ > 0.25. The higher order terms in 
the Ci are fixed to the SM values while the leading terms A±, with the exception of A 7 , are 
allowed to float. Since the branching fraction measurement of B — > X s j is consistent with the 



prediction within the SM, Aj is fixed at the SM value, —0.330, or the sign-flipped value, +0.330. 
We choose Ag/A-j and A10/A7 as fit parameters. The SM values for Ag and A10 are 4.069 and 
-4.213, respectively 9 . To extract the these ratios, we perform an unbinned maximum likelihood 
fit to the events in the signal box with a probability density function (PDF) that includes the 
normalized double differential decay width. 

We measure the integrated asymmetry .4fb> which is defined as 

- = f f \sgnjcos 8)g(q 2 ,9)d cos 9dq 2 

FB / l\g(q 2 ,o)d cos edg* 

We determine the yield in each q 2 and forward-backward regions from a fit to the Mb c distribu- 
tion. Then we correct the efficiency and obtain 

Afb(B -» K*£ + £~) = 0.50 ±0.15 ±0.02, 
Afb{B + -> K + £ + £~) = 0.10 ±0.14 ±0.01, (3) 

where the first error is statistical and the second is systematic. A large integrated asymmetry is 
observed for K*£ + £~ with a significance of 3.4c. The result for K + l + £~ is consistent with zero 
as expected. The fit results of ratios of Wilson coefficients are summarized in Table 1. Figure 2 
shows the fit results projected onto the background-subtracted forward-backward asymmetry 
distribution in bins of q 2 . 

Table 1: Ag/A-j and fit results for negative and positive A7 values. The first error is statistical and the 



second is systematic. 
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Figure 2: Fit result for the negative At solution (solid) projected onto the background subtracted forward- 
backward asymmetry, and forward-backward asymmetry curves for several input parameters, including the effects 
of efficiency; A7 positive case (A7 = 0.330, Ag = 4.069, A10 = —4.213) (dashed), A10 positive case (A7 = —0.280, 
A 9 = 2.419, A10 = 1.317) (dot-dashed) and both A 7 and A 10 positive case (A 7 = 0.280, A 9 = 2.219, A w = 3.817) 
(dotted). The new physics scenarios shown by the dot-dashed and dotted curves are excluded. 

The fit results are consistent with the SM values ^9/^7 = —12.3 and A\q/Aj = 12.8. In 
Fig. 3, we show confidence level (CL) contours in the (Ag/Aj, A\q/At) plane based on the fit 
likelihood smeared by the systematic error, which is assumed to have a Gaussian distribution. 
We also calculate an interval in AgAio/A 2 at the 95% CL for the allowed A 7 region, 

-14.0 x 10 2 < AgA^/A 2 - < -26.4. (4) 

From this, the sign of ^9^10 must be negative, and the solutions in quadrants I and III of Fig. 3 
are excluded at 98.2% confidence level. Since solutions in both quadrants II and IV are allowed, 



we cannot determine the sign of AjA\q. Figure 2 shows the comparison between the fit results 
for the negative A 7 value projected onto the forward-backward asymmetry, and the forward- 
backward asymmetry distributions for several input parameters. We exclude the new physics 
scenarios shown by the dotted and dot-dashed curves, which have a positive AgA\o value. 




Figure 3: Confidence level contours for negative Aj. Curves show ler to 5 a contours. The symbols show the fit 

(circle), SM (triangle), and ^4io-positive (star) cases. 



3 Evidence of the Purely Leptonic Decay B — > r v T 

We use a 414 fb _1 data sample containing 447 x 10 6 B meson pairs collected with the Belle 
detector. We use a detailed MC simulation, which fully describes the detector geometry and 
response based on GEANT 10 , to determine the signal selection efficiency and study the back- 
ground. The B~ — > t~v t signal decay is generated by the EvtGen package 11 . To model the 
background from e + e~ — > BB and continuum qq (q = u, d, s, c) production processes, large BB 
and qq MC samples corresponding to about twice the data sample are used. 

We fully reconstruct one of the B mesons in the event, referred to hereafter as the tag side 
(-E>tag), and compare properties of the remaining particle(s), referred to as the signal side (B s i g ), 
to those expected for signal and background. In the events where a B tag is reconstructed, we 
search for decays of B s \ g into a r and a neutrino. Candidate events are required to have one 
or three charged track(s) on the signal side with the total charge being opposite to that of 
-Btag- The r lepton is identified in the five decay modes, n'v^Vr, e~v e u T , ir~v T , Tt~it Q v T and 
ir~ n + n~ v T , which taken together correspond to 81% of all r decays 4 . The muon, electron and 
charged pion candidates are selected based on information from particle identification devices. 
For all modes except r~ — ► 7r _ 7r°^ r , we reject events with 7r° mesons on the signal side. 

The most powerful variable for separating signal and background is the remaining energy 
in the ECL, denoted as Eecl, which is sum of the energy of photons that are not associated 
with either the B tag or the 7r° candidate from the r _ — > tv~tv°i/ t decay. For signal events, Eecl 
must be either zero or a small value arising from beam background hits, therefore, signal events 
peak at low Eecl- On the other hand background events are distributed toward higher Eecl 
due to the contribution from additional neutral clusters. The Eecl signal region is optimized 
for each r decay mode based on the MC simulation, and is defined by Eecl < 0.2 GeV for the 
fi~i>nV T , e~u e v T and n~v T modes, and -Eecl < 0.3 GeV for the -k~it g v t and k~ tx~ v T modes. 
The -Eecl sideband region is defined by 0.4 GeV < -Eecl < 1-2 GeV for the vTv^Vt, e~v e v T 
and tt~v t modes, and by 0.45 GeV < Eecl < 1-2 GeV for the 7r _ 7r°z/ T and 7r~7r + 7r~^ T modes. 
Table 2 shows the number of events found in the sideband region for data {Nfg) and for the 
background MC simulation (N^^). Table 2 also shows the number of the background MC events 
in the signal region (-A^^p). In order to validate the Eecl simulation, we use a control sample of 
events (double tagged events), where the E tag is fully reconstructed as described above and B s { g 
is reconstructed in the decay chain, B~ — > D*°£~i> (D*° — > -D°7r°), followed by D° — > K~ir + or 



K 7r 7r + 7r + where £ is a muon or electron. Figure 4 shows the Eecl distribution in the control 
sample for data and the MC simulation scaled to equivalent integrated luminosity in data. 




Figure 4: -Ebcl distribution for the both B tagged events, where one B is fully reconstructed and the other B 
is reconstructed as B~ — > D*°l~v. The dots with errors indicate the data. The solid histogram represents the 
background from BB MC (B + B~ + B°B°), and the dashed histo gram shows the contribution from B°B° events. 

After finalizing the signal selection criteria, the signal region is examined. Figure 5 shows 
the obtained £"ecl distribution when all r decay modes are combined. One can see a significant 
excess of events in the E'ecl signal region below Eecl < 0.25 GeV. Table 2 shows the number 
of events observed in the signal region (iV b s ) for each r decay mode. 
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Figure 5: -Eecl distributions in the data after all selection requirements have been applied. The data and 
background MC samples are represented by the points with error bars and the solid histogram, respectively. The 
solid curve shows the result of the fit with the sum of the signal shape (dashed) and background shape (dotted). 

We deduce the final results by fitting the obtained Eecl distributions to the sum of the 
expected signal and background shapes. Probability density functions (PDFs) for the signal 
fs(EECh) and for the background /;,(-Eecl) are constructed for each r decay mode from the 
MC simulation. The signal PDF is modeled as the sum of a Gaussian function, centered at 
Eecl = 0, and an exponential function. The background PDF, as determined from the MC 
simulation, is parametrized by a second-order polynomial. The results are listed in Table 2. 
The number of signal events in the signal region deduced from the fit (N s ) is 21.2±|;^ when all 
r decay modes are combined. Table 2 also gives the number of background events in the signal 
region deduced from the fit (iVjj), which is consistent with the expectation from the background 
MC simulation {N™?). 

The branching fractions are calculated as B = N s /{2 ■ e ■ N B + B -) where N B + B - is the 
number of Y(4S) — ► B + B~ events, assumed to be half of the number of produced B meson 
pairs. The efficiency is defined as e = e tag x e sel , where e tag is the tag reconstruction efficiency 
for events with B~ — > t~u t decays on the signal side, and e sel is the event selection efficiency 
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Table 2: The number of observed events in data in the sideband region (N°^ c ), number of background MC events 
in the sideband region (N^^.) and the signal region (A^ c ), number of observed events in data in the signal region 
(N ohB ), number of signal (N B ) and background (N h ) in the signal region determined by the fit, signal selection 
efficiencies (e scl ), extracted branching fraction (23) for B~ — > t~v t . The listed errors are statistical only. 



listed in Table 2. When all r decay modes are combined we obtain a branching fraction of 
(1.06lo;2g) x 10~ 4 . The branching fraction for each r decay mode is consistent within error as 
shown in Table 2. 

Systematic errors for the measured branching fraction are associated with the uncertainties 
in the number of B + B~, signal yields and efficiencies. The total fractional uncertainty of the 
combined measurement is I^lo^j and we measure the branching fraction to be 

B(B~ -> t-P t ) = (1.06l°;|(stat)lg;l(syst)) x 10" 4 . 

The significance is 4.0a when all r decay modes are combined, where the significance is defined as 
S = ^-21n(£ /£ m ax), where 

£ m ax an d jCq denote the maximum likelihood value and likelihood 
value obtained assuming zero signal events, respectively. Here the likelihood function from the 
fit is convolved with a Gaussian systematic error function in order to include the systematic 
uncertainty in the signal yield. 
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